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Abstract 

Context. Apart from the discovery of hundreds of exoplanets, the high-precision photometry from the GoRoT and Kepler 
satellites has led to measurements of surface rotation periods for tens of thousands of stars, which can potentially be 
used to infer stellar ages via gyrochronology. 

Aims. Our main goal is to derive ages of thousands of field stars using consistent rotation period measurements derived 
by different methods. Multiple rotation periods are interpreted as surface differential rotation (DR). We study the 
dependence of DR with rotation period and effective temperature. 

Methods. We re-analyze the sample of 24,124 Kepler stars from Reinhold et al. (2013) using different approaches based 
on the Lomb-Scargle periodogram. Each quarter (Q1-Q14) is treated individually using a prewhitening approach. 
Additionally, the full time series, and different segments thereof are analyzed. 

Results. For more than 18,500 stars our results are consistent with the rotation periods from McQuillan et al. (2014). 
Thereof, more than 12,300 stars show multiple significant peaks, which we interpret as DR. Dependencies of the DR 
with rotation period and effective temperature, as shown in Reinhold et al. (2013), could be confirmed, e.g. the relative 
DR increases with rotation period. Gyrochronology ages between 100 Myr and 10 Gyr were derived for more than 17,000 
stars using different gyrochronology relations, most of them with uncertainties dominated by period variations. We find 
a bimodal age distribution for Teff between 3200-4700 K. The derived ages reveal an empirical activity-age relation 
using photometric variability as stellar activity proxy. Additionally, we found 1079 stars with extremely stable (mostly 
short) periods. Half of these periods may be associated with rotation stabilized by non-eclipsing companions, the other 
half might be due to pulsations. 

Conclusions. The derived gyrochronology ages are well constrained since more than ~93.0% of the stars seem to be 
younger than the Sun where calibration is most reliable. Explaining the bimodality in the age distribution is challenging, 
and limits accurate stellar age predictions. The relation between activity and age is interesting, and requires further 
investigation. The existence of cool stars with almost constant rotation period over more than three years of observation 
might be explained by synchronization with stellar companions, or a dynamo mechanism keeping the spot configurations 
extremely stable. 
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1. Introduction 

Encouraged by high-precision photometry of the CoRoT 
and Kepler satellites, measurements of stellar rotation pe¬ 
riods have become numerous in recent years, unprecedented 


et al. 2012| |McQuillan et al. 2013a| 

Nielsen et al. 

12013 

McQuillan et al.|2013bt |Walkowicz & ' 

3asri|2013 Reinhold| 


lAffer 


et al. 2013 McQuillan et al. 2014). Owing to magnetic brak- 
ing, the rotation period of a star is linked to its age. Stellar 
winds carry away charged particles along the magnetic field 
lines. As time proceeds, more and more material is dis¬ 
sipated, making the star slower and slower. This angular 
momentum loss over the stars’ lifetime results in the fact 
that young stars rotate faster than old ones, on average. 

A relation b etween stella r age a nd rotation period was 


first shown by Skumanichl (1972). This author further 


demonstrated that the level of chromospheric activity de- 
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pends on the stellar age. Succeeding measurements at the 
Mount Wilson Observatory corroborated the relations be¬ 
tween chromospheric activity index the rotation pe- 


riod, and the stellar age ( Noyes et al.jl984 Soderblom et al. 
1991). Additionally, since the 198Qs fast rotators are known 
to exhibit enhanced X-ray activity ( Pallavicini et al^|1981 


Pizzolato et al.| 2003). These coherencies are now estab¬ 
lished as activity-rotation-age relation. 

The stellar age cannot be measured directly but has to 
be inferred from other quantities. In principle, activity-age 
relations can be used to infer stellar ages, although their 
accuracy suffers from secular changes of the activity level. A 
classical method to infer stellar ages is isochrone modeling. 
The ages of stellar clusters can b e inferred from this method 
(see, e.g., Perryman et al.||1998 ), provided some knowledge 
of effective temperature, luminosity, and metallicity. This 
method can also be applied to field stars, however with 
typical uncertainties on the order of 20-30%. 

Over the past years, gyrochronology has become a 
promising method to derive stellar ages from their rotation 


1 













































Timo Reinhold, Laurent Gizon: Rotation, differential rotation, and gyrochronology of active Kepler stars 


periods (Barnes||2QQ3 2007 Mamajek fc Hillenbrand||2008 


Meibom et al. 2009 

Collier Cameron et al. |2009t James 

et al. 

|2010t Barnes 

|2010| IDelorme et al. |2011| Meibom 


et al. 


2011b). Barnes (2003112007) collected rotation peri- 


ods from open clusters of different ages, and identified two 
sequences in the color-period diagram, namely the /- and 
C-sequence. Using known cluster ages, this author derived 
empirical fits to the color (mass) dependence of the data 
on the /-sequence, according to 




( 1 ) 


The function g(t) oc is the rotation-age relation from 
Skumanich (1972), and f{B — V) an empirical fit to the 
{B — V)o color dependence of the rotation period Pj of the 
/-sequence stars refiecting their braking efficiency. {B — V)o 
colors are used as a substitute for stellar mass instead. 
Unfor tunately, gyrochronolo gy is poorly calibrated for old 


stars, 
the 1 


Meibom et al. 


(2011a) measured rotation periods for 
Gyr old cluster NGC6811 in the Kepler field. The 
oldest star used for calibration w as the Sun with an ag e of 


4.55 Gyr. Bridging this large gap Meibom et al. (2015) re¬ 
cently measured rotation periods of 30 stars in the 2.5 Gyr 
old cluster NGC6819, deriving a well-defined period-age- 
mass relation. At later ages wide bi naries might help con- 
straining the age-rotation relation (Chaname & Ramirez 

Another method to dete rmine stellar ages is p rovided by 


asteroseismology (see, e.g., Chaplin et al. 2014). Recently, 


efforts are underway to calibrate asteroseismology ages by 


comparing them to gyrochronology measurements (Garcia 
et al.||2014| |do Nascimento et al. 2014[ [Lebreton fc GoupII 


2014| [Angus et al.||2015D. Furthermore, [Vidotto et al. 
( |2014 ) found a correlation between the average large-scale 
magnetic field strength and the stellar age according to 
{\By\) oc similar to Skumanich’s law, calling this 

method magnetochronology. A re view article on var ious age 
dating methods was provided by Soderblom (2010). 

Gyrochronology relies on age calibration using open 
clusters of different ages, assuming their stars to be co¬ 
eval. Calibration becomes difficult for clusters older than 
^ 1 Gyr, though. One reason is that rotation periods are 
difficult to measure for slowly rotating stars. Furthermore, 
clusters become disrupted at that age (the degree strongly 
depends on the initial stellar cluster density), render¬ 
ing the cluster membership of individual stars uncertain. 
Additionally, gyrochronology relations are calibrated only 
for main sequence dwarfs. Bulk rotation period meas ure¬ 
ments might be largely cont aminated by subgiants (van 
Saders fc Pinsonneanlt 2Q13[ ), which have started evolv¬ 
ing off the main sequence. Thus, applying gyrochronology 
relations using subgiant rotation periods mi ght lead to a 
mis-c lassification of their evolutionary state (Dogan et al. 

2M^. 

Further uncertainties result from the fact that a star is 
by no means a rigid rotator with a well-defined rotation pe¬ 
riod. High-precision instruments like the CoRoT and Kepler 
telescopes provide the opportunity to observe multiple ro¬ 
tation periods associated with latitudinal differe ntial rota¬ 
tion, which was observed for many active stars (Reinhold 
et al.|2Q13 ). Moreover, spot rotation periods heavily rely on 
their evolutionary timescales, which becomes more impor¬ 
tant for less active stars. Both effects may contribute large 
uncertainties to the rotation period used in the color-period 
fits. 


Despite these drawbacks, gyrochronology relations sup¬ 
ply a straightforward way to infer stellar ages of field stars. 
Other approaches to measure stellar ages are usually ac¬ 
companied with large errors. Isochrone methods collapse 
for binary stars, which are assumed to be coeval, but do 
not appear on the same isochrone if the companions ex¬ 
hibit different masses. Asteroseismology ages strongly de¬ 
pend on some (usually unknown) model parameters (e.g., 
metallicity), which can lead to large errors. 

The knowledge of stellar ages is of utmost interest for 
galactic formation. We aim to provide stellar ages for thou¬ 
sands of stars, inferred from mean rotation period mea¬ 
surements via gyrochronology. Furthermore, we show that 
multiple significant periods are quite common among ac¬ 
tive stars, which we assign to surface differential rotation. 
These period fluctuations dominate the age uncertainties. 
As opposed to this, we also found a sub-sample of stars 
almost showing no period variations over more than three 
years of observation. 

The paper is organized as follows: Sect. presents the 
Kepler data products. The different approaches used to an¬ 
alyze the data are explained in Sect. Sect.JS contains 
the results, which are further discussed in Sect. We close 
with a brief summary of our results in Sect. 


2. Kepler data 


The Kepler satellite provides almost continuous observa¬ 
tions of the same field for more than four years (May 2, 
2009 - May 11, 2013). The data is delivered in quarters 
(Q0-Q17), each of ^ 90 d length (Q2-Q16), with exceptions 
for the commissioning phase QO (^ 10d), and Q1 (^33d). 
Unfortunately, observations stopped after one month of 
Q17 due to a failure of the third reaction wheel. This amaz¬ 
ing amount of Kepler data is publicly available and can be 
downloaded from the MAST archivqB 

Kepler data has been processed by different pipelines so 
far, starting with the Presearch Data Conditioning pipeline 
(PDC), designed to detect planetary signals. This pipeline 
was changed to th e so-called PDC-MAP (Maximum A 
Posteriori) pipeline ( Stumpe et al.||2012 Smith et ar]|20 12 ) 
because the PDC pipeline coarsely removed stellar variabil- 
ity signals. Recently, all Kepler data has been rep rocessed 
by the PDC -msMAP (multiscale MAP) pipeline (Stumpe 
et al. 2014). This new version applies a 20-days high-pass 
filter intending to detect smaller planets in the data. Thus, 
this pipeline version is not suitable to look for stellar vari¬ 
ability with a broad range of rotation periods since it di¬ 
minishes stellar signals of slow rotators. Unfortunately, data 
reduced by previous pipeline versions are not publicly avail¬ 
able a nymore. _ 

In Reinhold et al. (2013) we analyzed ^40,000 active 
sta rs with a variability ra nge i?var > 0.3% (for definition 
see Basri et al.|2010 2011), only using Q3 data. In ~ 24,000 
stars a clear rotation period was detected, and in ~ 18,000 
stars a second period was detected, which was assigned to 
surface differential rotation. Starting from this sample, we 
extend our analysis to all available data. Our goal is to 
detect consistent rotation periods throughout the quarters, 
and to refine previous differential rotation measurements, 
exploiting the much higher frequency resolution thanks to 
the longer time span. 


^ http://archive.stsci.edu/pub/kepler/lightcurves/tarfiles/ 
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We do not use the much larger sam ple of 34,030 stars 
with measured rotation periods from |McQuillau et ah] 
(2014) because we are mostly interested in measuring DR. 
Although 20,009 stars of our sam ple are contained in the 
sample of McQuillan et al. (2014), the remaining stars ei- 
ther do not belong to the periodic sample of |McQuillan| 
et al. (2014), or have a smaller average variability range 
(i^var < 0.3%), rendering them unsuitable for DR mea¬ 
surements. 

In this work we only use Q1-Q14 PDC-MAP data (ver¬ 
sion 2.1 for Q1-Q4, Q9-Q11 and version 3.0 for Q5-Q8, 
Q12-Q14). QO data was discarded because of its short time 
span and the much lower number of monitored targets. 
From Q15 on, only PDC-msMAP data was available, which 
is unsuitable for our purposes as explained above. 

Since we are only interested in rotation-induced stellar 
variability, we have to exclude targets showing other kinds 
of periodic variability. To reduce the number of false posi¬ 
tives, i.e., periodic variability not related to stellar rotation, 
we discarded 17 eclipsing binarie^ 878 planetary candi- 


datej^ 2 R R Lyrae stars (Kolenberg et al. 2Q1Q| 


et al. |2Q1Q| [Benko et ^|2Q1Q[ |Gnggenberger et al. 


Szabo 


Nemec et al.|2Qll||Moskalik et al.|2012||Molnar et ah 2012 


2012 


Nemec et al. 2013D, and 84 7 Doradns and J Sent! stars 


([Tkachenko et arf2013| Ulnsoy et al. 2014 Balona et al. 


2011 jUytterhoeven et al. 2011[|Lampens et al.|2013[|Balona 


2014). In total, 981 stars were discarded, leaving 23,143 


targets, which are analyzed as described in the following 
section. 


3. Methods 

After excluding binarity stars and pulsators, we are inter¬ 
ested in the stability of the rotational modulation. Many 
effects can change the shape of a light curve dominated 
by star spots rotating in and out of view, e.g., spots 
are created, others disappear while rotation takes place. 
The number of spots and their sizes are usually unknown. 
Differential rotation further hampers the detection of stable 
rotation periods. Sun spots change their preferred latitude 
of occurrence during the solar activity cycle, and therefore 
their rotation rate. But cyclic variability is also expected 
in other active stars. Moreover, Kepler suffers from instru¬ 
mental effects, which are corrected by the pipeline automat¬ 
ically. Improper correction can mimic long-term (periodic) 
variability. The Kepler satellite rolls between consecutive 
quarters to re-orient ate its solar arrays. Hence, in each quar¬ 
ter stars fall on different CCDs with different sensitivity. 

To account for all these effects, we analyze the data of 
the same star in different ways. First, we apply our anal¬ 
ysis to each quarter individually. By comparing periods of 
individual quarters among themselves, we shrink our initial 
sample to find stable rotation periods among many quar¬ 
ters. After that, we stitch together data from all quarters 
for each star in our sample, and apply a slightly different 
period search to the full light curve. In a last step, we ana¬ 
lyze different segments of the full light curve, because some 
effects, (e.g., differential rotation) are more visible in cer¬ 
tain segments rather than in the full time series. By com¬ 
paring the periods returned by each method, we hope to 


^ http://keplerebs.villanova.edu/ 

^ http: / /archive.stsci. edu/kepler/koi/search.php 


detect stable rotation periods, which are unlikely caused 
by instrumental systematics. 


3.1. Analysis of individual quarters 


To detect periodic signals in the data we use the Lomb- 
Scargle p eriodogram in a prew hitening approach as de¬ 
scribed in 
individual 


Reinhold et al 


_ ^(2013). Each quarter is analyzed 

y and in the same way. The analysis method in 
briefly summarized below. For details we refer the reader 
to the paper mentioned above. 

For each star and each quarter we compute the variabil¬ 
ity range i?var,Q- If ^var,Q > 0.3% holds for a certain quar¬ 
ter, we compute five Lomb-Scargle periodograms in a suc¬ 
cessive prewhitening approach. Since computing the Lomb- 
Scargle periodogram is equivalent to fitting a sine wave to 
the data, each prewhitening step yields a set of sine fit pa¬ 
rameters (period P/^, amplitude a/c, phase 0 /^, and a total 
offset c) for k = 1,...,5. The parameters which belong to 
the highest peak in the periodogram yield the best param¬ 
eters to fit the data. We use the returned values as initial 
parameters for a global sine fit to the data (t, y) according 
to 


yf'^t ^ sm(^ t - 0 /c) + c, 

k=l ^ 


( 2 ) 


This global sine wave is fit to the data through x^- 
minimization. To save computation time each light curve 
was binned to two hours cadences. The period associated 
to the highest peak in the Lomb-Scargle periodogram is 
called Pi, and represents the most significant period in the 
data. 

In some cases spots are located on opposite sides of the 
star, and only half of the true rotation period is detected. 
To reduce this number of alias periods we compare Pi to 
the set of periods Pk and their corresponding peak heights 
hk := h{Pk). If it holds 

|A-P/c/2| < 0.05Pi and hk>0.5hi, (3) 

then the period Pk is likely the true rotation period. In 
case Pi already was the correct period, we check if there 
are periods Pk satisfying 

|A-2P/c| < 0.05Pi and hk>0.5hi. (4) 

If that is not the case, the period Pk is used as rotation 
period, and we call this rotation period Pi. 

Since our primary goal is to detect differential rotation, 
we look for periods adjacent to Pi. A probable second pe¬ 
riod should satisfy 

0.01 < |Pi -P/cl/Pi < 0.30 (5) 

The period Pk with the second highest power in the 
prewhitening process satisfying Eq. is called P 2 . In the 
following section we compare period measurements of indi¬ 
vidual quarters among themselves. 


3.1.1. Comparison of periods from different quarters 

As discussed at the beginning of Sect.|^ the measured pe¬ 
riods Pi may change from quarter to quarter. That is espe¬ 
cially true for stars exhibiting a second period P 2 because 
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Quarter 

Pi [dJ 

P 2 [dj 

JJvar [%J 

QOl 

5.51 

4.98 

0.69 

Q02 

5.71 

6.06 

1.41 

Q03 

5.39 

5.14 

1.55 

Q04 

5.47 

6.04 

2.34 

QOS 

5.34 

5.74 

1.78 

Q06 

5.48 

5.12 

1.63 

Q07 

5.44 

5.20 

1.86 

QOS 

5.30 

5.81 

1.78 

Q09 

5.38 

5.79 

1.44 

QIO 

5.33 

5.74 

1.32 

Qll 

5.37 

6.02 

1.79 

Q12 

5.29 

5.89 

1.66 

Q13 

5.37 

5.95 

1.66 

Q14 

5.64 

5.26 

1.17 


Table 1. Period measurements of the star KIC 1163579 
from Fig[2 For all quarters Q1-Q14 the star was monitored, 
and the variability range i?var,Q > 0.3 %. In each quarter a 
second period was found, and no alias period was detected. 
For the star KIC 1163579 we find (Pi) = 5.43d. 


spots may have changed their location, size, or occurrence 
at all. Besides that, the period P 2 might be a spurious de¬ 
tection in some quarters. As a first step, we look for stars 
with stable rotation periods throughout the quarters. To 
’’pick and choose” the best stars shrinks our primary sam¬ 
ple, but provides a more reliable measure of the mean stellar 
rotation period, and eventually the stars’ DR. 

For each star and each quarter Q = 1,..., 14 we compute 
the variability rang^Pvar,Q- If ^var,Q > 0.3% we apply 
our prewhitening approach, which yields a period Pi,q. If 
^var,Q < 0.3%, or the star was not monitored in a certain 
quarter, the period Pi^q is set to zero. We only consider 
periods satisfying 0.5 d < Pi,q < 45 d. The upper limit ac¬ 
counts for the fact that we want to see at least two full 
rotation cycles during the limited quarter length of ^ 90 
days. The lower limit should exclude pulsating stars. We 
compute the relative deviation of the periods Pi^g from 
their median according |Pi,q — Pi|/Pi. Fast rotators with 
a median period Pi < 10 d are allowed to differ by 10 %, 
stars with longer periods by 20%. From all periods Pi^q 
satisfying these criteria, we define a mean rotation period 
(Pi). If more than 75% of the periods Pi^q satisfy this 
criterion, the star belongs to the so-called ’’good” sample. 
The ’’good” stars are shown in Fig. [^compared to previous 
measurements for the whole sample. We also checked for 
alias periods among our set Pi^q according to Eq. but 
without comparing any peak heights. Periods identified as 
such were discarded. The star KIC 1163579, randomly cho¬ 
sen from our sample, is shown in Fig. and its briefing 
is shown in Table In the following section we concate¬ 
nate data from individual quarters to analyze the full light 
curve. 


periods longer than 45 days, and 2 ) to obtain a higher fre¬ 
quency resolution in the Lomb-Scargle periodogram, since 
the peak width scales with the inverse of the time span. 

The first point is difficult to address. As mentioned ear¬ 
lier long-term instrumental effects are sometimes difficult to 
distinguish from slow rotation, especially in an automated 
period search of a large sample of objects. Furthermore, for 
this particular sa mple rotation periods less than 45 days 
were measured in Reinhold et al. (2013). Measurements of 
alias periods or spurious detections are possible but ex¬ 
pected to be rare. Nevertheless, we use an upper period 
limit of 60 days for the analysis of the full time series. 

The second point is the more interesting one because 
a higher frequency resolution offers the detection of indi¬ 
vidual peaks lying much closer than in the periodograms 
of the individual quarters. Thus, the frequency resolution 
is crucial for the detection of small values of differential 
rotation! 

To achieve an almost continuous time series of each star 
we chose the easiest way to concatenate consecutive quar¬ 
ters by dividing the light curves of each quarter by its me¬ 
dian and su btracting unity . Data outliers were removed as 
described in Garcia et al. (2011). Light curves exhibiting 
low-frequency trends, i.e., increasing or decreasing behav¬ 
ior over the full 90 days window, were considered as not 
properly reduced by the PDC-MAP pipeline, and therefore 
discarded. 

The method we apply to the full time series slightly dif¬ 
fers from the previously used sine fit approach. We only 
compute one Lomb-Scargle periodogram of the whole time 
series. Due to the largely increased amount of data each 
light curve was binned to six hours cadences to save com¬ 
putation time. This sampling still enables us to detect short 
periods down to half a day according to the Nyquist fre¬ 
quency /Nyq ~ 2 d~^. We identify the twenty highest peaks 
and search for periods within the limits 0.5 d < Pi < 60 d. 
Furthermore, we force a lower peak height limit oihi >0.10 
to get some significance for Pi, and check for alias periods 
according to Eq. and 

An important point to make is that we do not over- 
sample the periodog ram, in contrast to previous work 
( Reinhold et al.|[2M3 ). Eor the analysis of individual quar- 
ters oversampling was necessary because many cases only 
revealed a single broadened peak, rather than two or more 
distinct peaks. Therefore, a fine frequency sampling was 
necessary to subtract the correct period in the prewhiten¬ 
ing to be able to detect more than one period. Owing to the 
much higher frequency resolution here, the periodogram is 
able to reveal individual peaks. 

In the following section we try to assign a significance to 
individual peaks. The goal is to find out which peaks really 
carry information about different spot rotation periods, and 
which are related to spurious detections and/or stochastic 
effects, e.g., emerging and waning active regions. 


3.2. Analysis of the full time series 

In the previous section we only considered periods shorter 
than 45 days, due to the limited quarter length of ~90 
days. Stitching together data from all available quarters is 
the only way to achieve two major goals: 1 ) the detection of 

^ In the following, the variability range Pvar is meant to be 
the median of the periods Pi,q for the quarters Q1-Q14. 


3.2.1. Identification of significant peaks 

Analogous to Eq.|5]we search for periods Pk within 30% of 
Pi. We sort the peak heights in descending order according 
to hi > h 2 > ... and compute the so-called peak height ra¬ 
tio (PHR) hk/hk^i. This method is illustrated based on the 
example star KIC 1163579 in Eig. The idea is to detect 
peaks with comparable heights, and to see where the peak 
height drops to a significantly lower value from one peak to 
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Figure 1. Light curves (black) and sine fits (red) according to Eq. |^of the star KIC 1163579 for all quarters Q1-Q14. 
Period measurements of this star are given in Table 


the other. Hence, we compute the median of the PHR and 
search for the maximum deviation from this value. The re¬ 
lated index k = kmax yields the number of significant peaks. 
In Fig. the median is shown as solid red line, and the 
dashed red lines indicate the Tier region. For this example, 
it is evident that for k^ax = 3 the PHR is at maximum, 
which means that the peak height of P 3 is much higher than 
that of P 4 , compared to all peak heights within Pi ± 30 %. 
All periods with k < /Cmax are considered as significant, 
the others are discarded. The lower panel of Fig. shows 
the periodogram of the active star KIC 1163579 with the 
significant peaks found by this method marked in red. 

3.3. Analysis of segments of the full time series 

This section contains another approach we used to extract 
information about different periods in the data. The pri¬ 
mary peak in the periodogram of the full time series repre¬ 
sents the best sine fit period over the full observing time. 
This period must be interpreted as an average rotation pe¬ 
riod, because it fits certain parts of the light curve better 
than others. Adjacent peaks with less power are also needed 
to properly fit the data, and these periods may be more vis¬ 
ible in certain segments than in the full light curve. 

Minor peaks adjacent to Pi (e.g., peaks marked by blue 
asterisks in the lower panel of Fig. 1^ result from various 
phenomena: differential rotation, spot evolution, and data 
reduction. Differentiating between these effects is almost 



Figures. The peak height ratio (PHR) hk/hk^i versus 
peak number k of the periodogram from Fig. The solid 
red line shows the median PHR, and the dashed red lines 
mark the Tier levels. At k = 3 the deviation from the me¬ 
dian is at maximum, meaning that three significant peaks 
were found by this method. These are marked by red aster¬ 
isks in Fig. 

impossible, but we assume that DR is the dominating ef¬ 
fect in most cases. Spot evolution might be the strongest 
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KIC 1163579 



200 400 600 800 1000 1200 1400 

BJD - 2454833 [d] 



500 550 600 650 700 

BJD - 2454833 [d] 



Figure 2. Upper panel: Full light curve of the star KIC 1163579 from quarters Q1-Q14. Middle panel: Zoom into the 
above light curve. Periodicity and the existence of a second period is clearly visible, revealed by the double-dip structure 
around 550 days. Lower panel: Periodogram of the full light curve. Significant peaks are marked in red, other peaks in 
the range of Pi ± 30 % (indicated by the dashed blue lines) are marked in blue. 


contributor of spurious detections of DR, where multiple 
peaks cannot be associated to spots rotating at different 
latitudes. Analyzing different segments thus becomes im¬ 
portant, since the evolution of individual spots does not 
occur continuously in time, but may affect certain parts of 
the full light curve stronger than others. Another problem 
might be caused by stitching together individual quarters, 
where the end of one quarter and the beginning of the con¬ 
secutive one do not agree with the overall variability pat¬ 
tern. This problem is intrinsic to the PDC-MAP pipeline, 
which was designed to remove instrumental effects from in¬ 
dividual quarters, and not to conserve the overall variability 
pattern over the total observing time. 

To overcome these effects we analyze segments with dif¬ 
ferent lengths of the full light curve, aiming to detect peri¬ 
ods that are stable over different time scales. In the follow¬ 
ing, segments are defined as concatenated quarters, starting 
with data from Q1-Q2, and adding the subsequent quar¬ 
ter in the next step, i.e., the second segment contains data 
from Q1-Q3, and so on, with the last segment being the 
full light curve Q1-Q14. We compute the periodogram of 
each segmenlj^ interpolate it onto the frequency grid of the 


periodogram of the full light curv^ and add up the log¬ 
arithmic powers of all segments. This leads to extremely 
clear periodograms, likely canceling spurious detections. 

Fig. 13 shows periodograms of the defined segments of 
the star KIC 1163579, clearly revealing multiple periods 
between 5-6 days in all segments. The main period crys¬ 
tallizes out around 5.4 days. Each periodogram reveals a 
period around 6.1 days, which was also found in the pe¬ 
riodogram of the full time series (s. Fig. [^, but was not 
considered as a significant period by our method. Thus, 
the measured value of DR might be largely underestimated 
in this case. To account for such periods with minor power, 
we sum up the logarithmic power of all segments, which 
is shown in the upper panel of Fig. The main power is 
visible around Pi = 5.4 days, and also the first and sec¬ 
ond harmonic being the half and third of this period. In 
certain cases it occurs that the first harmonic, Pi/2, has 
the second highest or even the highest power in the pe¬ 
riodogram. Since we are interested in periods adjacent to 
Pi we subtract a fourth order polynomial (solid blue line 
in Fig. 1^, and exponentiate the difference, which is shown 
in the Lower panel of Fig. This procedure leads to ex¬ 
tremely clear periodograms. Periods lying with ±30 % of 


^ For t he individual segments we use the normalization from 
Eq. 22 in Zechmeister & Kiirster (2009). 


This is necessary because periodograms of different segments 
possess a different frequency resolution. 
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Figure 4. Periodograms of the segments of KIC 1163579. Periodicity is clearly visible between 4.5-6.5 days. The location 
of the main peak around 5.4 days does not change much, whereas the number of minor peaks increases toward longer 
segments, owing to the higher frequency resolution. 


Pi with powers > 0.5 hi are considered as significant. The 
dashed red line indicates the arbitrarily chosen significance 
limit. The period around 6.1 days survives this procedure, 
and will contribute to the DR measure in Sect. 14.21 Results 
combining our three different approaches are presented in 
the following section. 


4. Results 


We present rotation periods of more than 18 ,500 stars using 
the different approaches from Sects. [3TH3.3| Bas ed on these 
periods we search for differential rotation in Sect. |4.2[ Mean 
rotation periods together with their uncertainties are used 
to infer stellar ages through different gyrochronology rela¬ 
tions in Sect. |4.3j Finally, Sect. |4.4] is dedicated to a small 
fraction of stars exhibiting rotation periods very stable in 
time. 




4.1. Rotation periods 

In Fig. we present the mean rotation period (Pi), av¬ 
eraged over the quarters Q1-Q14, versus color (P — 17)0. 
The b l ack do ts show previous measurements from |Reinhold 


et al. (2013), only using Q3 data. Green dots mark the 


18,691 so-called ’’good” stars (s. Sect. 3.1), and the 454 
red and 625 blue dots indicate very stable rotators, which 
are discussed separately in Sect. |4.4[ The dashed blue lines 


Figure 5. Upper panel: Summed up logarithmic powers of 
the periodograms of the individual segments from Fig. 
The solid blue curve shows a fourth order polynomial fit. 
Lower panel: Exponential power of the upper panel with 
the polynomial fit subtracted. Peaks above the arbitrarily 
chosen significance limit of 0.5 (dashed red line) are con¬ 
sidered as significant, and used as a measure of DR in the 
following. In this picture the harmonics (Pi/2 and Pi/3) 
are also more pronounced. 
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No, "good" stars: 18,691 
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Figure 6. Mea n rotation periods (P i) averaged over the quarters Q1-Q14 plotted against color {B — V)o- Previous 
measurements (Reinhold et al. only using Q3 data are shown in black, ’’good” measurements surviving the criteria 

described in Sect. m are s hoW in green. Red and blue dots show stars with periods ve ry stable in time, which are 
discussed separately in Sect. 4.4 The dashed blue lines show isochrones from Barnes (2007), and the blue star marks the 
position of the Sun. 


show isochrones from Barnes (2007) and the blue star marks 
the position of the Sun. Around {B — V)o = 0.4 rotational 
braking due to magnetized winds becomes efficient. The 
4500 Myr isochrone acts as an upper envelope to our mea¬ 
surements up to (B — K)o — 1.0. Coeval stars redder than 
{B — V)o = 1.0 with supposably longer periods are missing 
due to the limited quarter length of ^ 90 days. A lower 
envelope to the rotation period distribution is given by 
the 200 Myr isochrone. Stars below this curve are either 
younger than 200 Myr, or are not suitable for the use of 
gyrochronology. The latter is also true for stars bluer than 
(B - V)o = 0.4 (s. Sect. [IXT] ). A second group of stars 
immediately leaps to the eye with short periods between 
0.5-2 days and {B — K)o < 0.4. Most of these stars exhibit 
perio ds v ery stable in time, which we discuss separately in 
Sect. |4.4| Comparing the black and green dots it is evi¬ 
dent that the period measurements have been considerably 
improved by incorporating more data. 

For the three different approaches (Sects. |3.1 


3.3) we 


compared ou r measurements to the st ate-of-the-art rotation 
periods from McQuillan et al. (2014) in Fig. In general, 
all methods show very good agreement. In each panel the 
dashed red line shows the one-to-one period ratio, and the 
upper and lower dashed blue lines indicate the 2:1 and 1:2 
period ratios, respectively. The left panel shows average pe¬ 


riods from the quarters Q1-Q14, using the green, red, and 
blue dots from Fig.[^ We find 17,674 stars matching the two 
samples. Thereof, more than 98.9 % of our measurements lie 


within 10% of those from McQuillan et al. (2014). A small 


fraction of 0.5 % alias periods was found. The middle panel 
shows rotation periods derived from the analysis of the full 
light curve. 15,082 stars are matching, 97.4% thereof do 
not differ by more than 10%, and 1.4% alias periods were 
found. The right panel shows measurements from the indi¬ 
vidual segments. 16,443 stars are matching, thereof 96.0 % 
within 10%, and 2.3% alias periods were detected. 

Combined period measurements from the methods 
above are shown in Fig. We computed the median of 
the periods Pi derived from each method, allowing for a 
median absolute deviation (MAD) of one day for periods 
shorter than twenty days, and a MAD of two days for longer 
periods. I n total, we find 1 8,599 stars matching with the 
sample of McQuillan et al. (12014 ). 97.6 % of the periods lie 
within 10 % of each other. In contrast to McQuillan et al. 


(2014) our sample contains stars hotter than 6500 K, and 
also a few stars which do not belong to their periodic sam¬ 
ple. Since we are searching for DR in the following sec- 
tion, we do not restrict our sample to stars matching with 
McQuillan et al. (2014), but consider all stars with mea- 
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Figure 7. Comparison of rotation periods Pi to the results of McQuillan et al. (2014) using different methods. The 
dashed red line denotes the 1:1 ratio, and the upper and lower dashed blue lines indicate the 2:1 and 1:2 period ratios, 
respectively. For each method the number of stars matching all criteria is given in the lower right corner of each panel. 


Combined measurements 



Figure 8. Rotation period measurements combining the 
three different approaches from Fig.[7l and compar ing them 
to AutoACF periods from|McQuillan et al.| (|2Q14|). 


sured period Pi within the above MAD limits, additionally 
satisfying log^ >3.5 and i?var > 0.3%. 


4.2. Differential rotation 

As shown in the previous section, rotation periods can be 
detected in a straight forward way by picking the highest 
periodogram peak and applying certain selection criteria. 
The situation is different when searching for differential ro¬ 
tation, acting as a perturbation to the main rotation pe¬ 
riod. We usually interpret the detection of a second pe¬ 
riod in the periodogram analysis as an indication of DR. 
Unfortunately, each method from Sects. 3.1 3. 3| yields dif¬ 
ferent results. For a certain star it occurs that one method 
returns a second period, whereas the other one does not. 
Even if all three methods yield a second period for a certain 
star, the periods found can differ a lot, depending on the 
different selection criteria and the frequency resolution. For 
that reason we define mean values for the DR of a certain 
star in the following. 

We pick the minimum and maximum of the set of sig¬ 
nificant periods, individually for Sects. and |3.3[ If both 
methods yield a minimum and maximum period for a cer¬ 


tain star, we compute the mean value, and call these periods 
Pmin and Pmax 5 respectively. We refrain from calculating 
minimum and maximum periods from individual quarter 
measurements owing to their lower frequency resolution. 
Thus, we define the relative and absolute horizontal shear 
^ • (-^max -frnin)/-fmax and dU :— 27r (1/Pmin iZ-^max)? 
respectively, as a measure of DR. In the following, we show 
how these two quantities correlate with rotation period and 
effective temperature. We only consider stars exhibiting a 
variability range Pyar > 0-3 %, which w as defined as a lower 
activity limit in Reinhold et al. (2013). 

Fig.i shows the relative shear a as a function of the 
minimum period Pmin- For stars cooler than 6700 K we find 
that a increases with rotation period. A contrary behavior 
is found for hot stars (Teff > 6700 K) populating the upper 
left corner {a > 0.02 and Pmin < 2d). These short period 
stars spread a wide range of a, and clearly do not follow 
the overall trend. We suggest that multiple period measure¬ 
ments in these stars might actually be due to pulsations or 
rapid spot evolution, and should not be interpreted as DR. 
Owing to magnetic braking, cool stars exhibit longer pe¬ 
riods, on average, thus populating the upper right part of 
Fig. 1^ In the lower left corner {a < 0.01 and Pmin < 3d) 
a mixture of all temperatures is found, indicating young 
fast rotating stars. Again, we warn the reader that these 
small a values can also be mis-classified as spot evolution. 
The dashed blue area s hows theoretical predictions from 
Kiiker & Riidigerl (2011) for models with O. 3 M 0 (bottom) 
to I.IMq (top). These simulations agree well with our ob¬ 
servations of relative shear increasing with rotation period. 
Furthermore, the hot stars are not covered by the model 
predictions, supporting our conclusion above. Almost 78 % 
of the fast rotators (Pmin < 2 d) exhibit rotation periods 
very stable in time. These very stable periods are consid¬ 
ered separately in Sect. |4.4[ 

We compare d our findings to previous measurements 
from Hall| (1991) who also found an increase of the relative 
shear witn rotation period. The dash-dotted and dashed 
black lines show linear fits to our measurements in log-log 
space yielding a relation a oc P^^^^ with c = 0.71 discarding 
the hot stars (blue data points) and c = 0.55 using all data 
points, res pectively. T he first value is consistent with the 
result from [Hali ( |1991 ) who found c = 0.79 ± 0.06. 

We also compared our results to measurements from 


Donahue et al. (1996). These authors found a relation be- 
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Figure 9. Minimum period Pmin versus relative differen¬ 
tial rotation a. The colors represent different temperature 
bins, the dotted black lines marks our detection limits. 
The dashed blue area s hows theoretical predictions from 
Kiiker & Rudiger I (|20 111) for models with 0.3Mq (bottom) 
to 1.1 Mq (top). We ffna that a increases with rotation pe¬ 
riod for stars cooler than 6700 K. The black lines show lin¬ 
ear fits to our measurements, discarding stars hotter than 
6700 K (dash-dotted line) and using all data points (dashed 
line), respectively. 


tween the mean rotation period (P) and the observed pe¬ 
riod spread AP according to AP oc (p^i-3±o.i^ Keeping in 
mind that a oc AP/(P), we find AP oc (P)^+^ = (P)^-^^. 
T his value is slightly bi gger than the value of 1.3 found 
by Donahue et al.| ( 1996| ), regardless of whether discarding 
or keeping hot stars. The discrepancy might be explained 
by the fact that these authors only considered FGK stars, 
whereas we also incorporate M stars in our sample, usually 
exhibiting large values of a. 

Fig. shows the dependence of the relative shear 
on effective temperature for a subset of stars satisfying 
2d < Pi < 3d. The dash-dotted red and dashed orange 
curves show theoretical predictions from [Kiiker fc: Rndiger| 


revised temperatures from Huber et al. ([2014 ), being on 
average 200 K hotter for hot stars and 200 K cooler for cool 
stars, compared to previous temperatures from the KIC, 
which might explain this offset. We do not want to stress 
the comparison of our observations to the models because 
we do not find a clear trend when using al l rotation peri¬ 


ods. T his is consistent with previous work (Reinhold et al. 
2013), where only a shallow trend of a towards cooler stars 


has been found. We now interpret this shallow increase of a 
as a consequence of the increase of rotation period towards 
cooler stars. 

In Fig.[^we plot the absolute shear dfl against rotation 
period Pmin- We find that dft does not strongly depend on 
rotation period over a wide period range. Towards fast ro¬ 
tators with periods on the order of a few days, the absolute 
shear increases, although showing large scatter. Again, we 


(2011) using a model with an equatorial rotation period of 
2.5 days. Between 3000-6000 K the red curve matches the 
observations quite well, although the slope of the curve is 
too small. Stars hotter than 6000 K are well represented 
by the slope of the orange curve, although our observa¬ 
tions are offset towards hig her temperatu r es. W e are using 


Figure 10. Teff versus a for stars satisfying 2d < Pi < 3d. 
The dash-dotted red and dash ed orange curve show models 
from Kiiker & Rudiger (2011). The dotted black line marks 
the upper limit of a. No distinct correlation between Teff 
and a was found. 


find that the upper left corner is populated by hot stars (s. 
Fig.§, which are clearly separated from the overall trend. 
Again, we warn the reader that these large dft values might 
not be associated with strong surface shear, but with dif¬ 
ferent pulsation frequencies or rapid spot evolution. The 
dashed blue area s hows theoretical models from Fig. 3 in 
Kiiker & Riidiger (2011) with 0.3Mq (bottom) to I.IMq 
(top), which cover most of our measurements, but do not 
touch the hot stars with dft > 0.1 radd“^. 

Equivalent to Fig. we applied a linear fit to our mea¬ 
surements yielding dQ oc (P)^ oc with c = —0.29 

discarding hot stars (dash-dotted line) and c = —0.45 us¬ 
ing all data points (dashed line) in Fig.[^ res pectively. We 
find go od agreement with measnrementsirom [Barnes et al. | 
(2005) claiming dfl oc (11)confirming the weak de- 
pendence of the absolute shear on the rotation rate. 

The dependence of the absolute shear on effective tem¬ 
perature is shown in Fig.f^ Our measurements (gray dots) 
suggest the existence oiTwo distinct regions with differ¬ 
ent behavior of dO. From 3500-6000 K the absolute shear 
slightly increases showing weak dependence on tempera¬ 
ture. Above 6000 K dH steeply increases with temp erature, 
although sh owing large scatter. Data collected by [Barnes 


et al.[ (2005) is pl otted as purple diamond s, and the corre- 
spending fit from Collier Cameron[ ( [2007 ) is shown as dot¬ 
ted purple line. These authors predicted a very strong tem¬ 
perature dependence (dH oc T^ff^), which is not supported 
by our meas urements. Light blue diamonds sho w measure¬ 
ment^ from Ammler-von Eiff & Reiners (2012[) that are in 
good agreement with our fin dings. Theoretical predictions 
from [Kiiker fc Rudiger (2011) are shown as dash-dotted red 
and (lashed orange curve. These authors found no sufficient 
match when fitting their results with one polynomial over 
the whole temperature range. Therefore, they suggested a 
different behavior of dQ above and below ^6000K. Their 
curves fit remarkably well with the mean values of our mea¬ 
surements, shown as thick blue line. The matching gets 


^ The values were taken from Table 2 in [Ammler-von Eiff fc[ 
Reiners (2012) assuming an inclination of 90°. 


10 


































































Timo Reinhold, Laurent Gizon: Rotation, differential rotation, and gyrochronology of active Kepler stars 


No. stars: 12,319 
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Figure 12. Effe ctive temperature Te g versus absolute shear dQ su mmarizing different meas urements: Purple diamonds 
~ ' 10^, the dotted purple curve from|Collier Cameron| (|2007|). Light blue diamonds show 

m & Reiners I (|2012|). Our measurements are shown as gray dots. The thick blue line 


were taken from Barnes et a h 
measurements from I Ammler- 


von 


represents the weighted mean of our measure ments for 250 K temperatu re bins. The dash-dotted red line and the dashed 
orange line show theoretical predictions from Kiiker & Rudiger (2011). The dotted black line represents our detection 
limit, and the red star marks the position of the Sun. We find that dil shows weak dependence on effective temperature 
between 3300-6200 K. Above 6200 K dQ strongly increases showing large scatter. 


even better when using th e old KIC temperat ures, rather 
than the new values from Huber et al. 

200 K hotter (co oler) for th e hot 


__ (2014), which are 

co ol) st ars, on average. 


The model from |Kuker & Rudiger (2011) is drawn from 
a star rotating with an equatorial period of Peq = 2.5 d. 
Our observations contain various rotation periods, which 
might explain the spread. In Sect. we discuss reasons for 
the observed spread, including cases with multiple periods 
mis-interpreted as DR. The measured values of Pmin, Pmax, 
Of, and dft are collected in Table The measured period 
differences are used as uncertainties for determining gy¬ 
rochronology ages in the next section. 


4.3. Stellar Ages 
4.3.1. Gyrochronology 

After measuring mean rotation periods P and period varl 
ations AP for two-thirds of the sample, we calculate stel 
lar ages t along with their un certainties At using differ¬ 


ent gyrochronology relations ([Barnes _[2QQ7 Mamajek & 


KIC 

TeS 

(K) 

Rva.r 

(%) 

Pmin 

(d) 

■fmax 

(d) 

a 

dQ 

rad/d 

9771948 

5678 

1.63 

15.115 

15.689 

0.037 

0.015 

11122123 

5960 

1.17 

9.883 

10.955 

0.098 

0.062 

2558314 

4157 

0.82 

24.790 

28.171 

0.120 

0.030 

12317549 

5793 

1.06 

7.377 

7.996 

0.077 

0.066 

3103637 

5332 

0.71 

26.372 

33.500 

0.213 

0.051 
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The full 


Table 2. DR measurements used in Figs 
table is available in machine-readable form in the online 
journal. 


viated as B07, MH08, and MMS09, respectively. Ages are 
calculated according to Eq. 3 in Barnes (2007): 


log t = — [log P — log a — b log X ], 
n 


(6) 


with X := (B — Ejo — c and different fit parameters a^h^c^n 
(see Table 1^. Colors (P — E)o are obtained b y converting 
g — rpl to (P — V) color using the relation from Jester et al. 
(|2005f and subtracting the excess B — V reddening E{B — 


Hillenbrand 2008[ Meibom et al. 2009p , hereafter abbre- ® Taken from the Kepler Input Catalog (KIC). 
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Figure 11 . Minimum period Pmin versus absolute differen¬ 
tial rotation dfl. The colors represent different temperature 
bins, the dotted black lines mark our detection limits. The 
dashed blue area shows theoretical predictions from |Kuker| 
& Rudiger (2011). We found that dft is almost independent 
of rotation period over a wide period range, but slightly in¬ 
creases towards short periods. Stars hotter than 6700 K do 
not follow the overall trend. The black lines show linear fits 
to our measurements, discarding stars hotter than 6700 K 
(dash-dotted line) and using all data points (dashed line), 
respectively. 


At=- J{AP/PY + [hAX/Xy + de^, 

n 


(7) 


with AP := (Pmax ~ Pmin)/2 being the uncertainty of the 
mean period P, and AX b eing the {B — V) diff erenc e be¬ 
tween the co nversion from [Jester et al. ( 2005 ) and Bilir 

li 


et al 


Table (compare Eq. 10 in |Barnes]|2007 ). If DR was 
3 ted the period uncertainty dominates the age uncer- 


(2005|). Frequently, this difference is rather small 
and no uncertainties for the SDSS g — r colors were found 
in the literature. Thus, we set a minimum uncertainty of 
±0.01 mag to the {B — V) colors with smaller uncertainties. 
The term contains uncertain ties of the fit parameters 
from 

detected the period uncertainty 
tainty, especially for the slowly rotating stars with AP on 
the order of a few days. 

Distributions of the derived ages are shown in Fig. 

For our stellar ages sample we retain stars that are not 
contained in the McQuillan et al. (2014) san iple , but remove 
those with stable rotation perioas (s. Sect. |4.4| ). Th e latter 
migh t not have spun down to the /-sequence yet (s. [Barnes 
2003 for terminology), or their rotation might be controlled 
by non-eclipsing companions. Additionally, we tighten our 
limit of the surface gravity to log^ > 4.2 to ensure that 
gyrochronology relations are only applied to dwarf stars. 


Ref. 

a 

b 

c 

n 

o 

1 

p[dj ■ 

B07 

0.7725 

0.601 

0.400 

0.5189 

> 0.42 

> 1.5 ‘ 

MHOS 

0.4070 

0.325 

0.495 

0.5660 

> 0.51 

> 1.5 ^ 

MMS09 

0.7700 

0.553 

0.472 

0.5200 

> 0.50 

> 1.5 < 


Table 3. Fit parameters and period and color constraints 
used in the different gyrochronology relations. 


V )^. Age uncert ainties are calculated according to Eq. 10 
in [BarnS l( |2007t : 


Figure 13. Stellar age distri butions derive d from differ¬ 
ent gyrochronology relations: |Barnes| (|2QQ 7|); [Mamajek & 
Hillenbrand (2008); [Meibom et al. ( 2Q09p . Tn eir calibra- 
tion ranges are given in Table [^ Stars with derived ages 
younger than 100 Myr (left of the dashed black line) should 
be treated with caution. 


The MHOS and MMS09 distributions contain less stars 
than the B07 distribution owing to the different color ranges 
(see Table 1^. The color and period constraints are needed 
to ensure that our field star sample obeys the same (or 
at least a similar) dependence of rotation period on color 
-and age as the cluster stars used for calibration. The de¬ 
rived ages are collected in Table [^ and their reliability is 
discussed in more detail in Sect. [51 

Derived ages younger than 100 Myr should be treated 
with caution. These stars cover the full {B — V)o range, but 
mostly exhibit rotation periods less than five days. They 
likely have not yet converged to the /-sequence, and thus 
are not suitable for applying gyrochronology relations. The 
right edge of the distribution with derived ages older than 
10 Gyr can also not be trusted since roughly half of the stars 
exhibit ages older than the universe. These values result 
from a combination of long rotation periods (P > 20) days 
and effective temperatures > 5500 K, with either or both 
of them being erroneous. 

Using the B07 distribution 17,623 stars possess ages 
between 100 Myr and 10 Gyr. Ther eof, 90.7% are you nger 
than 4 Gyr, in good agreement with Matt et al. (|2015|) esti¬ 
matin g ^ 95 % comparing the sample from McQuillan et al. 
(2014) to model predictions. Less than 0.62% of the de- 
rived ages are greater than 10 Gyr, and less than 2.2% of 
the B07 stars lie in the critical calibration region younger 
than 100 Myr, providing some confidence in the derived age 
distribution. 


4.3.2. Activity-age relation 


Inspired by Fig. 8 in 
interested in deriving 


Soderblom et al. 
similar 


( [199ip we are 
relation. 


a similar activity-age 
Unfortunately, spectra of Kepler stars are lacking, so we 
cannot compare the derived ages to the established chromo¬ 
spheric activity measure Phk* Nevertheless, the variability 
range Pvar can be used as an activity indicator in a sta¬ 
tistical sense. Using the B07 distribution we plot the age 
against the variability range in Fig. [T^ which was inspired 


12 













































































































Timo Reinhold, Laurent Gizon: Rotation, differential rotation, and gyrochronology of active Kepler stars 


KIG 

-Rvar 

(%) 

TeS 

(K) 

{B-V)o 

(mag) 

A{B-V)o 

(mag) 

p 

(d) 

AP 

(d) 

tBOT 

(Myr) 

AiB07 

(Myr) 

iMH08 

(Myr) 

AiMHOS 

(Myr) 

iM09 

(Myr) 

AiM09 

(Myr) 

flag 

7799260 

2.14 

5437 

0.671 

0.010 

23.837 

0.000 

- 

- 

- 

- 

- 

- 

ev 

11135449 

1.20 

4731 

1.032 

0.010 

31.783 

0.775 

2197 

162 

3154 

364 

2370 

313 

- 

3751292 

4.21 

6089 

0.556 

0.010 

3.656 

0.000 

172 

16 

241 

35 

279 

48 

- 

8481446 

0.39 

5699 

0.628 

0.010 

15.115 

0.000 

1710 

134 

1892 

222 

2213 

259 

- 

3757951 

6.17 

5345 

0.677 

0.010 

1.918 

0.006 

- 

- 

- 

- 

- 

- 

vs 


Table 4. Rotation periods, — colors, and gyrochronology ages using different calibrations. If no ages are provided, 
either the period or the color lies outside the valid calibration range (see Table Super-stable and very stable stars 
(flagged as ”ss” and ”vs”, respectively) were discarded from the stellar ages sample. Furthermore, stars with log^ < 4.2 
may have evolved off the main sequence (flagged as ”ev”). For these stars we provide a rotation period but no ages since 
gyrochronology relations are not calibrated for evolved stars. Further table entries are the KIC number, the variability 
range i?var, and the effective temperature Teff. 




Figure 14. Gyrochronology ages plotted against the variability range for 500 K temperature bins between 3200-6200 K. 
The age distribution shows a bimodality between 3200-4700 K, which vanishes for hotter stars. The variability range 
i?var decreases with age, becoming more emphasized towards hotter stars (increasing slope of dashed red line). 


by Fig. 4 in McQuillan et al. (2014). The ages were de¬ 


rived onl y using periods lying c 
found by 
lower rig' 


_ oser than 10% to periods 

found by |McQnillan et ah (2014). From the upper left to the 
ight, the temperature increases from 3200-6200 K in 
500 K intervals. The upper panels (3200-4700 K) shows a 
bimodality of the age distribution, which vanis hes for hot¬ 
ter st a rs. Th is incident was first detected by |McQuill^ 


et al. (2013a) for the Kepler M dwarf periods, and con- 
firmed later that the bimodality extends to hotter stars 
(Reinhold et al. 2013 McQuillan et al.|2014 ). Moreover, the 
gap separating the two peaks shifts towards younger stars 
with increasing temperature, starting at ~ 800 Myr for the 
coolest stars (3200-3700 K), and descends to 500 Myr for 
stars between 4200-4700 K. For each temperature bin the 
point distribution shows that the variability range decreases 
with age, although each distribution exhibits large scatter 
in i?var- This general behavior is expected from the obser¬ 
vation that young stars are, on average, more active than 


old ones. In each panel two point clouds are visible. Thus, 
we separated the two clouds into a so-called young and old 
sample. To guide the eye we empirically drew a dashed red 
line between the two samples, with increasing slope towards 
hotter temperatures, and varying offset in each panel. Stars 
lying below (above) the dashed red line belong to the young 
{old) sample, respectively. There are other ways of defin¬ 
ing a young and an old sample, e.g., using the horizontal 
dashed black line separating the two peaks of the bimodal 
distribution, but we wanted to emphasize the correlation 
between activity and age. We were curious if the parame¬ 
ters of these distinct samples substantially differ. Thus, we 
plotted histograms of common stellar parameters such as 
log^, FeH, brightness, and so on. Unfortunately, we found 
no major difference in the parameters that could explain 
the evidence of the two point clouds. However, we found a 
correlation between the ages shown and the corresponding 
peak heights of the primary periodogram peaks. The peak 
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heights of the young sample were, on average, higher than 
the peaks of the old sample. We interpret this result in 
terms of young stars being more active. Hence, their light 
curves are more sinusoidal, because they suffer less by DR, 
spot evolution, or instrumental flaws, all effects disrupting 
the light curve shape. 


4.4. Extremely stable periods 

From the analysis of the individual quarters we found pe¬ 
riods Pi,Q that are extremely stable in time. These stars 
have previously been shown in Fig. and are analyzed 
in more detail here. To quantify their temporal stability 
we computed the median absolute deviation MAD (Pi) := 
|Pi,Q ~ (Pi)!? categorized two groups of stable periodic 
stars: super-stable stars with MAD (Pi) < 0.001 d, and very 
stable stars satisfying 0.001 < MAD(Pi) < 0.01 d. These 
stars are flagged in the last column of Table We empha¬ 
size that these periods are stable over more than three years 
of observation! 

Figure shows the periods and effective temperatures 
of both groups. Super-stable and very stable periods are 
shown as red and blue dots, respectively. Inner green dots 
denote stars where a second period was found. The tem¬ 
perature and period distributions are shown in the upper 
and right panel, respectively. The period distribution shows 
that most of the super-stable stars exhibit short periods less 
than one day, with a mean period of (Pi) = 0.95 d. The very 
stable stars extend to longer periods up to ~ 12 d, with a 
mean of (Pi) = 1.90d. Stars with a second period were 
found in both groups, but only between 0.5-4.2 days. The 
temperature distribution shows that the vast majority of 
both groups exhibit temperatures less than 8000 K, with a 
mean of ^ 6800 K for the super-stable stars. The very sta¬ 
ble stars are, on average, TOOK cooler. Both distributions 
populate the full temperature range, but stars with a sec¬ 
ond period were almost exclusively found below 8000 K. 
Between 8000-10,000 K a dearth of second period stars was 
found, and only few second period stars were found above 
10,000 K. 

There exist two general explanations for the observed 
stable periodic variability. One possible explanation are 
stellar pulsations, which are known to be very stable in 
time. The other process serving as an astronomical clock 
is synchronization by a companion. We calculated the 
MAD (Pi) with the intention of disentangling these two 
processes. In the period regime of 0.5-4 days 6 Scuti, 
7 Dor, and hybrid pulsators thereof are expected. The 
boundaries of the so-called instability strip (i.e., the re¬ 
gion in the Hertzsprung-Russell diagram populated by 
these pulsators) is expected roughly between 6500-8800 K. 
Surprisingly, most of the stable stars exhibit temperatures 
less than 8000 K. A fraction of them might actually be 7 Dor 
stars between 6500-7000 K, but there is no pulsation mech¬ 
anism able to produce such stable periods in the cool stars 
regime. Stable periods above 8000 K are likely caused by 
pulsations because spots are not necessarily expected for 
such hot stars. Thus, we favor the conclusion that the pe¬ 
riodicity of stars cooler than 6500 K is caused by spots on 
the stellar surface stabilized by non-eclipsing companions, 
either due to interactions with another star or a close-in 
planet. Moreover, this hypothesis is supported by the ob¬ 


servation that stars with multiple periods (indicative for 
DR) were mostly found below 8000 K. 


5. Discussion 


5.1. Rotation 

— We measured rotation periods for a statistically mean¬ 
ingful ensemble of stars. In total, more than 18,500 ro¬ 
tation periods we derived using different approaches, 
all revealing very good agreement with the results from 


McQuillan et al. (2014). 

— As discussed in the previous section, we found 1079 pe¬ 
riods extremely stable in time, with a median absolute 
deviation less than 0.01 d. For stars cooler than 6500 K 
binarity is the favored explanation. In our total sample 
(see Table 0 we have 5124 stars with Pi < 10 d and 
Teff < 6500rK. 573 of these stars exhibit stable peri¬ 
ods, which corresponds to a per centage of ~ 11.2 % 


van 


Saders & Pinsonneault (2013) state that tidally syn- 
chronized binaries are fast rotators with periods less 
than 10 days and that they contaminat e the field with 
4%, compar ed to 11% in the Hyades (Duquennoy & 


Mayor 1 1991 ). Our rate is almost three tinies higher than 


the expected percentage of tidally locked field star bi¬ 
naries, but comparable to the expected percentage in 
the Hyades. One possible explanation may be the ex¬ 
istence of non-transiting planets, which are able to de¬ 
crease the a ngular momentum lo ss rate due to magne¬ 
tized winds (Cohen et al. 2010 ). Additionally, pulsat¬ 
ing stars cooler than 6500 K might contribute to this 
rate. Other expla nations f or sta ble rotation might be 
different dynamos. Brown (2014) suggests the so-called 
Metastable Dynamo^ where stars are born rapidly ro¬ 
tating with weak coupling to the wind. Other dynamo 
mechanisms might generate strong magnetic fields lead¬ 
ing to long spot lifetimes. 


5.2. Differential rotation 


— Exact values of a and dQ are hard to determine, and 
depend on the particular threshold used. Depending on 
which periods are selected as Pmin and Pmax, the abso¬ 
lute values of a and dO can differ a lot. Nevertheless, 
each method from Sects. 3.1f 3.3| produces the same cor¬ 
relations for a and dD with temperature and period. 
Combined measurements from the different approaches, 
as described at the beginning of Sect. |4.2[ were used to 
provide average values of Pmin and Pmax- As discussed 
by the example of Fig.[^the observed spread in dQ can 
only partially be exp lained by theory. |Collier Camer^ 


& Donati ( 2002 ) and Donati et al. (2003) found a large 
spread in ranging from 0.046 — 0.091 rad d“^ for the 
active star AB Dor. Evolving spot configurations might 
be an explanation. Hence, all measurements here should 
be interpreted in a statistical way. 

— Although the derived values are method dependent, the 
statistical averages are not. Our results are in good 


agreement with pr evious measurements ( Hall|_ iMT] 
Donahue et al.|1 99^ and theoretical predictions ^Kiiker 
&: Rudiger||201ip. 


— Simulations of spotted stars exhibiting rapid spot evolu¬ 
tion are able to generate beat-shaped light curves, mul¬ 
tiple periodogram peaks, and therefore able to mimic 
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Figure 15. Stars with super-stable periods 
(MAD(Pi) < 0.001 d, red), and very sta¬ 
ble periods (0.001 < MAD(Pi) < 0.01 d, 
blue). The inner green dots denote stars 
where a second period was found. The 
temperature and period distributions are 
shown in the top and right panel, respec¬ 
tively. 


DR. This interpretation was not considered so far. We 
think that spot evolution may play a role, but we have 
no way to discriminate between these two phenomena. 


5.3. Stellar Ages 

— All gyrochronology relations have been calibrated by 
ground-based observations of open clusters and Mount 
Wilson stars. The stars used in the different calibrations 
exhibit different period and color ranges. Age calibra¬ 
tion was performed using the Sun as an age anchor. 
Thus, the relations are not tested for stars older than 
the Sun. Applying these relations to stars with a wider 
period and color range might lead to less accurate ages, 
mean which 

— Gyrochronology ages are most reliable between 500- 
2500 Myr. Depending on their braking efficiency some 
stars younger than 500 Myr may not have converged to 
the /-sequence yet. Thus, their periods may not be suit¬ 
able for the use of gyrochronology. The calibration be¬ 
comes even worse for stars with derived ages less than 
100 Myr. Such stars might belong to the (7-sequence, 
obeying a physically different behavior. We do not trust 
derived ages younger than 100 Myr or older than 10 Gyr. 

— Subgiants and main seq uence stars obey a different 
rotation-age relationship (Garcia et al. 2014). In this 
study we attempt to exclude evolved stars by setting 
a lower limit to the surface gravity of log^ > 4.2. 
Unfortunately, Kepler does not provide stellar lumi¬ 
nosity classification, so our sample might be contami¬ 
nated by subgiants. Contamination mig ht be as large 
as 35% for field sta rs as pointed out by van Saders fc 
Pinsonneault| (|201^. 

tn( 


range in the sense that young stars are more active, on 
average. Interestingly, this trend becomes more distinct 
towards hotter stars, as indicated by the dashed red line 
in Fig. which lacks an explanation so far. 

A comparison of gyrochronology and asteroseismol- 
ogy ages is ch allenging. Much pr ogress was recently 
made (see, e.g., Angus et al. 2015). Most rotation pe¬ 
riod and asteroseismology samples do not overlap be¬ 
cause strong activity, essential for achi eving spot ro¬ 
tation periods, damps mode excitation (Chaplin et al. 


2011). Furthermore, Kepler lacks bright stars, which 


are needed for asteroseismology. Upcoming missions will 
hopefully change this unpleasant situation in the near 
future. 

6. Summary 


We re - analyz ed the sample of 24,124 stars from Reinhold 


et al. 


( |2013| ) using Q1-Q14 data. Good agreement was 
foun d with previous rotation periods and measurements 
from McQuillan et al. (2014). We searched for deviations 


from the mean rotation period using the Lomb-Scargle peri- 
odogram in different approaches, aiming to detect multiple 
significant periods, and assigning them to surface differen¬ 
tial rotation. The general tr ends of the DR with period 
and temperature, observed in Reinhold et al. (2013), could 
all be confirmed although individual measurements of a 
and dU may differ due to the different frequency resolu¬ 
tion of the full time series and the 90-days time base of 
a single quarter. In general, the new measurements are in 


very good agreement with previous observations ( Hall[ 
Donahue et al .||1996) and theoretical predictions (Kuker & 
Rudiger|201ip . Stellar ages were derived from gyrochronol 


— The existence of tne two point clouds in Fig. for 
stars with 3200 < Teg < 4700 K is a matter of debate. 
McQuillan et al.| (2014) suggested that the bimodality 
of the age distribution can be understood in terms of 
two distinct star formation events in the solar neigh¬ 
borhood. Stellar ages are correlated with the variability 


ogy relations provided by different authors, with uncertain¬ 
ties that are dominated by the period spread. A bimodal 
age distribution was found between 3200-4700 K, vanish¬ 
ing for hotter star. The derived ages show a correlation 
with the variability range serving as an activity indicator. 
Furthermore, we found 1079 stars exhibiting a very stable 
period, with a median absolute deviation less than 0.01 d. 
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The almost constant periods of the hot stars may be ex¬ 
plained by pulsations, whereas the stability of the cooler 
star (Teff < 6500 K) may be explained with synchroniza¬ 
tion of the orbital period of a non-eclipsing companion. 
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